Optical wave turbulence and the condensation of light 
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In an optical experiment, we report a wave turbulence regime that, starting with weakly nonlinear 

waves with randomized phases, shows an inverse cascade of photons towards the lowest wavenumbers. 
We show that the cascade is induced by a six-wave resonant interaction process and is characterized 
by increasing nonlincarity. At low wavenumbers the nonlincarity becomes strong and leads to 
modulational instability developing into solitons, whose number is decreasing further along the 
beam. 

PACS numbers: 190.0190, 190.4223, 190.3270, 190.3100, 190.6135. 



The idea to create the state of optical wave turbu- 
lence (OWT) has excited the minds of scientists for over 
thirty years, and it was a subject of a rather large num- 
ber of theoretical papers [1, 2, 3, 4, 5]. Indeed there 
axe some far reaching fluid analogies in the dynamics of 
nonlinear light, for example vortex-like solutions [6, 7], 
shock waves [8] and weakly interacting random waves 
whose dynamics and statistics has similarities to the sys- 
tem of random waves on water surface [9, 10]. OWT 
was theoretically predicted to exhibit dual cascade prop- 
erties similar to 2D fluid turbulence, namely the energy 
cascading directly, from low to high frequencies, and the 
photons cascading inversely, toward the low energy states 
[1, 2, 3, 4, 10]. When the nonlinearity is small, OWT can 
be described by the theory of weak turbulence (WT) [10] 
which possesses classical attributes of the general turbu- 
lence theory, particularly predictions of the Kolmogorov- 
like cascade states, which in the WT context are called 
Kolmogorov-Zakharov (KZ) spectra. It appears that 
OWT has two KZ states: one describing the direct energy 
cascade from large to small scales, and the second one - 
an inverse cascade of wave action toward larger scales. 
The inverse cascade is particularly interesting because 
in the optics context it means condensation of photons. 
Furthermore, it was theoretically predicted, and numeri- 
cally observed in some WT systems, that in the course of 
the inverse cascade the nonlinearity will grow, which will 
eventually lead to the breakdown of the WT description 
at low wavenumbers and to the formation of coherent 
structures [1, 3, 5, 11, 12, 13]. In optics, these can be 
solitons or collapses for focusing nonlinearity or a quasi- 
uniform condensate and vortices for the defocusing case. 

The final thermalized state was studied extensively 
theoretically in various settings for non-integrable Hamil- 

tonian systems starting with the pioneering paper by Za- 
kharov et al [14], see also [15, 16, 17, 18, 19, 20, 21, 22, 
23, 24]. The final state with a single soliton and small- 
scale noise was interpreted as a statistical attractor, and 
an analogy was pointed out to the over-saturated vapor 



system, where the solitons are similar to droplets and 
the random waves are like molecules [24]. There, small 
droplets evaporate while the big ones gain size from the 
free molecules, resulting in the decrease in the number of 
droplets. In our work we put emphasis on turbulence, i.e. 
on a transient non-equilibrium process leading to ther- 
malization rather than the thermal equilibrium itself. 

We report the first experimental evidence of OWT: 
starting with weakly nonlinear incoherent waves we 
observe an inverse cascade of the photons to lower 
wavenumbers and growth of the nonlinearity leading to 
the formation of multiple coherent solitons, which further 
merge into a single strong soliton. This corresponds to 
Bose-Einstein condensation (BEG) of photons, which is 
an optical analogue of atomic BEG reported in [25, 26]. 

The motion of photons to different energy levels (and 

toward the lowest one corresponding to condensation) 
can be achieved via nonlinear interactions of photons, 
provided for example by the Kerr effect. The problem, 
however, is that the nonlinearity is usually very weak 
and it is an experimental challenge to make it overpower 
the dissipation. In our view this was the main obstacle 
with the photon condensation setup in a 2D Fabry-Perot 
cavity theoretically suggested in [27]. 

The key feature in our setup is that we traded one 
spatial dimension for a "time axis" , namely we consid- 
ered a time independent 2D light field where the prin- 
cipal direction of the light propagation plays the role of 
"time" . This allowed us to use a nematic liquid crystal, 
which provides a high level and tunable optical nonlin- 
earity [28, 29]. The slow relaxation time of the reorienta- 
tional dynamics is not a restriction in our setup because 
our system is steady in time. Similar experiments were 
first reported in [30], where it was shown that a beam 
propagating inside a nematic layer undergoes a strong 
self-focusing effect followed by filamentation, soliton for- 
mation, when the light intensity is increased. More re- 
cently, a renewed interest in the same setup has led to 
further studies on the optical solitons and modulational 
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instability regimes [31, 32, 33]. However, all these ex- 
periments have in common a large value of the input in- 
tensity, of the order of lO'^ W/cm'^, therefore the optical 
nonlinearity is rather high and, as soon as the beam enter 
the liquid crystal, the soliton or modulational instability 
regimes appear immediately, bypassing the WT regime. 

In our experiment, we fix a much lower input inten- 
sity, of the order of 10~^ W/cm'^, and, for the first time, 
we show that, before the modulational instability and the 
soliton formation, corresponding to the breakdown of the 
weak nonlinearity, there is a WT regime characterized 
by an inverse cascade of the photon number. Besides 
the weak nonlinearity, to allow for the wave-mixing and 
the inverse cascade development, we also need to prepare 
a proper initial condition: we start with weakly nonlin- 
ear incoherent waves, characterized by high wavenumbers 
and random phases. 

The ID system is rather different from the 2D and 
3D systems considered theoretically before, and in this 
paper, for the first time, we present a theory of ID 
OWT, which involves coexisting random waves (inter- 
acting via a six-wave process) and solitons (to which 
the random waves condense). We complement our study 
with the numerical simulation of the dynamical equation, 
and we observe that the theory, numerics and experiment 
are consistent and complementary with each other, giv- 
ing a rather complete description of the ID OWT phe- 
nomenon. 
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FIG. 1: (color online) . Experimental setup: a laminar shaped 
input beam propagates inside the liquid crystal (LC) layer; 
random space modulations are imposed at the entrance of 
the cell by means of a spatial light modulator (SLM). 



The experimental apparatus is shown in Fig. 1. It con- 
sists of a laminar shaped beam propagating inside a ne- 
matic liquid crystal (LC) layer. Before entering the cell, 
the beam passes through a spatial light modulator (SLM) 
that, through suitable intensity masks, modulates the in- 
put profile for injecting random phased fields with large 
wavenumbers. The LC cell is made by sandwiching a ne- 
matic layer (E48) of thickness d = 50 fj,m between two 
20 X 30 mrrL^ glass windows. On the interior, the glass 
walls are coated with Indium-Tin-Oxide (ITO) transpar- 
ent electrodes. We have pre-treated the ITO surfaces in 



order to align all the molecules parallel to the confining 
walls (nematic director along z). The LC layer behaves 
as a positive uniaxial medium, with rig = 1.7 the extraor- 
dinary and Ho = 1.5, the ordinary refractive indices [29]. 
LC molecules tend to turn more along the applied field 
and the refractive index n{9) follows the distribution of 
the tilt angle. We apply a 1 kHz electric field with rms 
voltage Vq — 2.5 V so that the molecular director is pre- 
set to an average tilt 6. When a linearly polarized beam is 
injected into the cell, the LC molecules reorient towards 
the direction of the input polarization. The input light 
comes from a diode pumped solid state laser, A = 473 
nm, linearly polarized along y and shaped as a thin lam- 
inar Gaussian beam of 30 iim thickness, 1.8 mm width. 
The total power was 200 iiW, corresponding to an input 
intensity of 370 mW/cm^, which is low enough to ensure 
the weakly nonlinear regime. The beam evolution inside 
the cell is monitored with an optical microscope and a 
CCD camera. 

Theoretically, the beam evolution is described by a 
propagation equation for the input beam coupled to a 
relaxation equation for the LC dynamics 
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2«9ir- + ir^ + kln^aip = 0, 
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where ip{x, z) is the complex amplitude of the input 
beam propagating along "time axis" z, x the coordi- 
nate across the beam, a{x^ z) the liquid crystal reori- 
entation angle, Ua = rie — Uo the birefringence of the 
LC, fco the optical wavenumber, the vacuum permit- 
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tivity and = {■nK/2ls.e) {d/Vo) the electrical coher- 
ence length of the LC [34], with K the elastic constant. 



— /cq (n^ + n^/2) and Ae the dielectric anisotropy. 



Note that fixes the typical dissipation scale, limiting 
the extent of the inertial range in which the WT cas- 
cade develops. In other contexts, see e.g. [31, 32, 33], 
such a spatial diffusion of the molecular deformation has 
been denoted as a nonlocal effect. In our experiment, for 
Vq — 2.5 V we have = 9 fim. The experimental setup 
corresponds to the long- wave regime, kl^ <JC 1, where k 
is the typical wavenumber inside the range of the inverse 
cascade. Therefore, in this paper we consider the limit 
kl(^ ^ 1 of Eqs. (1) and (2), which enables us to form a 
single equation for the input beam [35] : 
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Eq. (3) conserves the energy 
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and the number of photons, N = J \ip\'^ dx. In the 
limit £on^l'^kQ/2K 0, Eq. (3) becomes the Un- 
ear Schrodinger equation and has linear wave solutions 
ip{x,z) ^ flfe cxp(— ia)feZ + ifc.x) with "frequencies" lo = k"^ 
and constant complex amplitudes a^- For weak nonlin- 
earity the amplitude become weakly dependent on 
"time" z. Note that the leading nonlinear term cannot 
generate a cascade in k because it corresponds to an in- 
tegrable ID Nonlinear Schrodinger equation. Thus, we 
retain the sub-leading term and apply the WT theory, 
which allows us to describe weakly nonlinear waves with 
random phases of Ofc. To apply the WT theory, we also 
need to verify that the linear dynamics dominate in the 
system. The ratio of the linear term and the leading 
nonlinear term of Eq. (3) is 
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where / is the input intensity. For / = 370 mW^/cm^ 
we have J ~ 100, which ensures the validity of the WT 
approach. 

For the wave action spectrum Uk = (akal) (the averag- 
ing is over the random phases), the WT approach yields 
to the following kinetic equation [35] , 

-g^ = A / /fcl2345 S{k + fci + fc2 - fcs - fc4 - fcs) 

x5{u)k + coi + u)2 — — UI4 — LJ5)dk-idk2dksdk4 
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(wave action in this case) is needed to fill the infinite iner- 
tial range. (Indeed, the integral of Uk ~ k~^/^ converges 
at fc = 0). In these cases the turbulent systems have a 
long transient (on its way to the final thermal equilib- 
rium state) in which the scaling is of the KZ type. This 
is because the initial condition serves as a huge reservoir 
of the cascading invariant. Note that the situation here 
is not specific for WT only and it is valid generally for 
turbulence. For example, it is valid for Navier-Stokes tur- 
bulence, i.e. the Kolmogorov-Obukhov spectrum, which 
is also finite capacity. 

To set up the inverse cascade, we have injected pho- 
tons at small spatial scales with random phases. This is 
made by creating through the SLM a random distribution 
of diffusing spots with the average size 35 fim. The 
numerical initial condition is more idealized and strictly 
localized at a small-scale range: we excite five wavenum- 
bers with constant amplitude around ~ 1.5 x 10^ 
ram~^ , with the phase of ■0^ being random and inde- 
pendent at each k. Moreover, we apply a Gaussian fil- 
ter in physical space to achieve a beam profile compara- 
ble to that of the experiment. Applicability of the WT 
approach for the numerical simulation is verified by the 
calculation of parameter J, which agrees with the exper- 
iment and is of the order of J ~ 100. Experimentally, 
we measure the light intensity I{x, z) = I'i/'p and not the 
phases of ij) and, therefore, the spectrum rifc is not di- 
rectly accessible. Instead, we measure the spectrum of 
intensity N{k,z) = \Ik{z)\'^. The scaling for in the 
inverse cascade state is easy to obtain from the wave ac- 
tion spectrum (8) and the random phase condition; this 
gives 



and A = 9nk^nl^lfef)/2048K^ 

The kinetic equation has important exact power law 
solutions rift = C\k\~'' where C and i' are constants. In 
particular, u = and v = 2 correspond to thermody- 
namic equilibria with equipartitions of the particle den- 
sity and the energy respectively. In addition, there are 
KZ power law solutions corresponding to a direct energy 
cascade from low to high fc's, 

nk = C\k\-\ (7) 
and an inverse wave action cascade from high to low fc's. 
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Note, that a kinetic equation can be derived directly 
from Eqs. (1) and (2) without the assumption of kl^ <C 1. 
We will consider the general case, including the short- 
wave limit of Eqs. (1) and (2) in an extended publication 
[35]. 

Here, we concentrate on the inverse cascade. The in- 
verse cascade spectrum is of a finite capacity type, in a 
sense that only a finite amount of the cascading invariant 



Experimental and numerical spectra of the light inten- 
sity are shown in Figs. 2 and 3 respectively. In both cases 
one can see an inverse cascade excitation of the lower k 
states, and good agreement with the WT prediction of 
the intensity spectrum (9). 

Closeness of Eq. (3) to integrability means that we 
should expect not only random waves but also soliton- 
like coherent structures. In the inverse cascade setup the 
sohtons appear naturally. Indeed, the WT description 
(Eq. (6)) breaks down when the inverse cascade reaches 
some low k's [35]. Modulational instability develops at 
these scales, which results in filamentation of light and 
its condensation into coherent structures - solitons. 

Verification of the long-wave limit kl^ <C 1 and devia- 
tion from integrability is checked by considering the ratio 
of the two nonlinear terms in Eq. (3) , which is estimated 
in Fourier space as i? ~ k^l^. We find from Figs. 2 and 
3 that the inverse cascade is approximately in the region 
k ^ 10* — 10^ giving an estimation of i? ^ 10~^ — 1. 
Note, that if R is too small, then we are close to a purely 
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FIG. 4: Experimental results for intensity distribution I{x, z). 
Area marked by the dashed line is shown at a higher resolution 
(using a larger magnification objectif). 



FIG. 2: (color online). Experimental spectrum of the light 
intensity, Nt ~ |/fep at two different distances z. 
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FIG. 3: (color online). Numerical spectrum of the light in- 
tensity, Nk ~ |/fcp at two different distances z. Averaging is 
taken over a small finite time window and over ten realiza- 
tions. 



integrable system, which would be dominated by soh- 
tons and lack cascade dynamics, and if too high then our 
model breaks down. 

Example zooms of the intensity distribution I{x, z) 
showing the beam evolution during propagation in the 
experiment and in the numerics are displayed in Fig. 4 
and Fig. 5 respectively. In the high resolution insets on 
Figs. 4 and 5 one can visually observe that the typical 
scale increases along the beam which corresponds to an 
inverse cascade process. Furthermore, in both Fig. 4 and 
Fig. 5 one can see the formation of coherent solitons out 
of the random initial field, with the overall number of 
solitons reducing as the beam propagates. Experimen- 
tally, the condensation into one dominant soliton is well 
revealed by the linear intensity profiles I{x) taken at 
different propagation distances, as shown in Fig. 6 for 
z = 0.3, 4.5 and 7.5 mm. Note that the amplitude of the 
final dominant soliton is three orders of magnitude larger 
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FIG. 5: Numerical results for intensity distribution I{x,z). 
The frame on the left is a magnified section of the initial 
propagation of the beam. 



than the amplitude of the initial periodic modulation. 
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FIG. 6: (color online) Linear intensity profiles I{x) taken at 
different propagation distances, z = 0.3, 4.5 and 7.5 mm. 

We perform two numerical simulations, one with a do- 
main comparable with the experiment (Figs. 5 and 3) 
and another with a lower intensity initial condition over a 
longer propagation distance z, to allow us to clearly dis- 
tinguish between the WT inverse cascade and MI (Fig. 7) . 
Deviations from integrability result in the interactions of 
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FIG. 7: (color online). The k — uj spectrum of the wave field 
at 2 = 2.1 m. ujt ^ l/256qZ| and fc. = l/\/l28Ze. The 
Bogoliubov dispersion relation is shown by the solid line. 



the solitons with each other and with the random waves, 
which lead to changes of the soHton strengths. Both 
the experimental and the numerical results in Fig. 4 and 
Fig. 5 indicate that solitons gradually merge, reducing 
the total number. The observed increase of the scale and 
formation of coherent structures represent the condensa- 
tion of light. 

Separating the random wave and the coherent soliton 
components can be done via performing an additional 
Fourier transform with respect to "time" z over a finite 
z- window determined by the frequency scale [5] . Such nu- 
merically obtained (A;,w)-plot is shown in Fig. 7. There, 
the incoherent wave component is distributed around the 
wave dispersion relation, which is Bogoliubov-modified 
by the condensate [35] and is shown by a solid line in 
Fig. 7 (the line ends at a finite k below which uj be- 
comes imaginary corresponding to modulational insta- 
bility) . This distribution is narrow for large k which cor- 
responds to weak nonlinearity, and it gets wider toward 
low k, which coincides with the growth of the nonlinearity 
and the breakdown of the WT applicability conditions. 
For these low k values one can see pieces of slanted lines 
(under the dispersion curve) each corresponding to a co- 
herent soHton, whose speed is equal to the slope. 

In conclusion, we have presented the first experimental 
implementation, numerical simulations and theory of ID 
OWT. We observe an inverse cascade of photons toward 
the states with lower frequencies with a respective KZ 
spectrum predicted by the WT theory. The inverse cas- 
cade is accompanied by increasing nonlinearity and the 
eventual breakdown of WT, leading to the light condens- 
ing into coherent structures - solitons. 
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